Gaussian localizable entanglement 
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We investigate localization of entanglement of multimode Gaussian states into a pair of modes by 
local Gaussian measurements on the remaining modes and classical communication. We find that 
for pure states and for mixed symmetric states maximum entanglement between two modes can be 
localized by local homodyne detections, i.e. projections onto infinitely squeezed states. We also 
show that non-Gaussian measurements allow to localize more entanglement than Gaussian ones. 



PACS numbers: 03.67.-a, 42.50.-p 
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Quantum entanglement, the non-classical correlations 
exhibited by quantum systems, lies at the heart of quan- 
tum information theory. Of particular interest are multi- 
partite entangled states which provide a resource for one- 
way quantum computing [l[ and could form a backbone 
of the quantum communication network, where a part 
of the entangled system is located at each node of the 
network. For quantum communication purposes, it is of 
great interest to determine how much entanglement can 
be localized on average between two nodes of the network 
by performing local measurements on the states located 
at the remaining nodes of the network and announcing 
the measurement results to the two nodes. This so-called 
localizable entanglement has been introduced and stud- 
ied in the context of quantum spin chains 0, H[ ■ 

In the present paper we investigate localization of en- 
tanglement in quantum continuous- variable systems 
More specifically, we consider multimode Gaussian states 
and investigate how much entanglement can be localized 
between two modes by local Gaussian measurements on 
the other modes. We prove that for pure Gaussian states 
as well as for mixed symmetric Gaussian states the opti- 
mal strategy is to carry out a balanced homodyne detec- 
tion (BHD) on each mode, i.e. to project it on infinitely 
squeezed state. Interestingly, we find that if we allow 
for non-Gaussian measurement strategies then we can 
localize more entanglement. The entanglement localiza- 
tion can be demonstrated experimentally with present- 
day technology. Various multimode entangled Gaussian 
states of light can be generated by combining single- mode 
squeezed states on an array of beam splitters [f| and 
highly efficient homodyne detectors are also available. 

The entanglement properties of multimode Gaussian 
states have been investigated previously [1, @] with 
particular focus on the symmetric Gaussian states in- 
variant under arbitrary permutation of modes [j| [t| 
ITol . [TT| . These latter states could be used to estab- 
lish a continuous- variable quantum teleportation network 
0, [l2| where quantum teleportation occurs between two 
(arbitrarily chosen) modes A and B and the other par- 
ties assist the teleportation by performing local measure- 
ments Mj on their modes Cj and sending the outcomes 
to the receiver B. Adesso and Illuminati 0] determined 
the optimal multimode symmetric state that for a given 
total amount of squeezing maximizes fidelity of telepor- 



tation of coherent states from A to B and also loga- 
rithmic negativity of the effective two-mode state of A 
and B when each Mj is balanced homodyne detection 
of p quadrature. Here we prove for all mixed symmet- 
ric Gaussian states that such local balanced homodyn- 
ing is optimal among all (not only local) Gaussian mea- 
surements and maximizes the entanglement established 
between modes A and B by Gaussian measurements on 
modes Cj. 

Consider TV-mode Gaussian state pabc shared among 
N parties A, B, Cj, j = 1, . . . , N~2, with each party pos- 
sessing a single mode. Parties Cj attempt to increase the 
entanglement between A and B by making local Gaus- 
sian measurements and communicating the measurement 
outcomes to A and B. By a Gaussian measurement on 
mode Cj we mean any measurement consisting of using 
auxiliary modes prepared in vacuum states, passive and 
active linear optics (beam splitters, phase shifters and 
squeezers) and BHD. Any such measurement can be de- 
scribed by the positive operator valued measure (POVM) 



Uj(aj) = ^Djiaj^DUaj). 



(1) 



Here Dj(aj) — exp(<x,cj — ot*Cj) denotes the displace- 
ment operator, n° is a density matrix of a (generally 
mixed) single-mode Gaussian state with covariance ma- 
trix and zero displacement and a.j is a certain linear 
combination of the quadrature values measured by the 
BHDs. In particular, homodyne detection on Cj is re- 
covered in the limit of infinitely squeezed state n°. A 
crucial feature of the Gaussian measurement is that 7^ 
depends only on the structure of the linear optical net- 
work but not on the measurement outcomes of the BHDs. 
The normalization Tr[H°] = 1 implies that 



iy D j (a j )I^D}(a j )d 2 aj = l j , 



(2) 



which ensures the completness of the POVM ((T|). With- 
out loss of any generality we can assume that each n° 
is a projector onto a pure Gaussian state because any 
mixed Gaussian state can be expressed as a mixture of 
pure Gaussian states and classical mixing cannot increase 
entanglement. 

The elements of the total POVM describing measure- 
ment on all modes Cj can be written as a product of 



the single-site elements, Ilc(a) = (S)jLi Hj( a j)> where 
a = (ai, • • • , oln— 2)- The measurement outcome a is 
obtained with probability density P(a) = Tt[1ab <8> 
IIc(Q!)pasc] an d the resulting normalized density ma- 
trix of the conditional bipartite state shared by A and B 
reads 



cab (<*) = 



1 



P(a 



■Trc[&AB ®n c (a) pabc]- (3) 
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FIG. 1: Decomposition of pure three-mode Gaussian state. 



Let -E[ovib] denotes a measure of entanglement of a bi- 
partite state (Jab ■ We define the Gaussian localizable en- 
tanglement El,g between A and B as the maximum en- 
tanglement that can be, on average, established between 
two parties A and B by local Gaussian measurements 
performed by Cj on their modes and by communicating 
the measurement outcomes to A and B. We have 



E L , G = max 



P{a)E[a AB {a)]da, 



(4) 



where the maximum is taken over all local Gaussian 

povms n c . 

It follows from the properties of Gaussian operations 
and measurements jl 31 ] that for all measurement out- 
comes a the conditionally prepared state <jab{o<-) is a 
Gaussian state with fixed covariance matrix and varying 
displacement which depends linearly on a. The entan- 
glement properties of a Gaussian state depend only on 
the covariance matrix, because the displacement can be 
set to zero by means of suitable local displacement op- 
erations Dj. It holds that E[aAB( a )] — E[<jab(0)], Va 
and from the definition ^ it immediately follows that 
El g — maxn c E[ctab(0)]- In order to determine the 
Gaussian localizable entanglement it thus suffices to op- 
timize over all local projections onto pure single-mode 
squeezed vacuum states. Generally, this optimization is 
still a daunting task and can be performed only numer- 
ically. However, we shall see that in case of three-mode 
pure states a fully analytical expression can be derived. 
Moreover, for general multimode pure Gaussian states we 
show that the optimal measurement involves homodyne 
detection on each mode so it suffices to optimize over the 
phases 9j which specify the quadratures being measured. 

Consider a pure three- mode Gaussian state \"4>)abc 
shared by parties A, B and C and characterized by a 
covariance matrix jabc- After projection of mode C 
onto a pure Gaussian state, the modes A and B will be 
in a pure Gaussian state \4>) ab with covariance matrix 



1AB = 



7A Sab 
6 ab IB 



(5) 



where ja (ib) denotes covariance matrix of mode A (B) 
and Sab contains correlations between the quadratures 
of the two modes. A unique measure of entanglement of 
pure states is provided by the entropy of entanglement, 
which is the von Neumann entropy of the reduced density 
matrix of one party .EU^ab] = — Tr^ log 2 p A ], where 
Pa = Ttb[|0)(0|ab]- The von Neumann entropy of a 



single-mode Gaussian state with covariance matrix ja is 
a function of the symplectic invariant tia — ( ^/det 7 a — 
l)/2. Explicitly, the entropy of entanglement reads 



E\ 



1 AB\ 



(n A + l)log 2 (n A + 1) - n A log 2 (n A ). (6) 



In order to maximize -E[|</>)ab] we have to maximize 
det7A. 

We prove that it is optimal to measure a suitably cho- 



sen quadrature of mode C xc,e = (ce 10 



Consider a bipartite splitting AB\C of the pure state 
\'4 ) )abc- It has been shown [3] that there exist uni- 
tary Gaussian transformations U\ B and IA C acting on 
modes AB and C, respectively, which transform the 
three-mode pure Gaussian state \4>)abc into a product 
of a two-mode squeezed vacuum state in modes A and C 
and a vacuum state in mode B, hi 
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u. 



I ABC 



\\)ac\0 )b- H ere \X) A c = vT^E^o A» A |n)c 
A 2 = (Vdet7c - l)/(Vdet7c + 1), 7c is CM of mode C 
prior to measurement and \n) denotes the n-photon Fock 
state. The situation is illustrated in Fig. [1] The transfor- 
mation Vic can be absorbed into the Gaussian measure- 
ment on mode C. Projection of one part of the two-mode 
squeezed vacuum state onto single-mode squeezed state 
with covariance matrix 7c = W (0)V (r)W T (9) , where 



W(6) = 



cos ft 
— sin( 



sin f 

cos ( 



V(r) 









,-2r 



prepares the mode A in a similar pure single-mode 
squeezed vacuum state |s; 9) with covariance matrix ja = 
W{-9)V{s)W T {~9), where 



.,2s 



1-A 2 + (1 + A 2 )e 
1 + A 2 + (1-A 2 )e 2 



2r 



(8) 



The resulting two-mode Gaussian state of modes A and 
B can be obtained from the product Gaussian state 
of modes A and B by the action of Uab, \4>)ab — 
Uab\s',9)a\0)b- Let Sab denotes the symplectic matrix 
corresponding to the unitary Uab which in the Heisen- 
berg picture governs the linear transformation of quadra- 
ture operators. We decompose Sab with respect to the 
A\B splitting, 



Sai 



Saa Tab 
Tba Sbb 



(9) 



The covariance matrix ^ab of the state \4>) ab can be ex- 
pressed as jab = Sab{ia © Is) Sa b , where the identity 



matrix represents the covariance matrix of vacuum 
state. After a straightforward calculation we arrive at 
j a = Saa^aS^a + TabTab- To calculate det 7^ which 
we want to maximize we use the formula for determinant 
of a sum of two 2x2 symmetric matrices X and Y, 



det(X + Y) = dct X + dct Y + Tr[XRYR 7 



(10) 



where R 



1 

-1 



Taking into account that deter- 



minant of a covariance matrix of a pure Gaussian state 
is equal to unity, det 7,4 = 1, we arrive at 

det 7A = (det Saa) 2 + {detT AB ) 2 + Tr[%M], (11) 

where M — Saa-^abTab^ 7 Saa is a symmetric posi- 
tive semidefinite matrix. M can be diagonalized by or- 
thogonal rotation, W(9 )MW T (9 ) = di&g{M xx ,M pp ). 
Since we optimize over all phases 9 we can, without loss 
of any generality, make the substitution 9 i— ► 9 — 9q and 
assume that M is diagonal. The nontrivial part of det ^a 
which should be maximized then reads 

e 2s (M xx cos 2 6+Mpp sin 2 9)+e~ 2a (M xx sin 2 9+M pp cos 2 9) 
The maximum of this function should be found in the in- 



where e 2s " 



(1+A 2 )/(1-A 2 ), 



terval s G 

c.f. Eq. ([8|). The maximum squeezing ±s max is obtained 
if mode C is projected onto the eigenstate of quadra- 
ture operator using homodyne detection. It is straight- 
forward to check that for any 9 expression (TT21) is max- 
imized if s — s max or s — — s max . This proves that the 
optimal Gaussian measurement on C which maximizes 
the entanglement between A and B is homodyne detec- 
tion. This remains valid even if we take into account the 
squeezing operation Lie which is included in this mea- 
surement. The only difference is that instead of quadra- 
ture qc we should measure quadrature xc,9 opt defined 

as UcqcUc = a x x c + a pPc = \Ja 2 x + a 2 , x c ,e a 



The 



optimal phase O pt specifying the quadrature xc,9 opt can 
be determined analytically by solving quadratic equation 
for tan 9 op t [IB] ■ Alternatively, the three-mode pure state 
can be transformed by local unitary Gaussian operations 
to a standard form, where all correlations between am- 
plitude quadratures Xj and phase quadratures pk vanish 
[1 01 ] . In this case it is optimal to measure either xq or 
Pc- 

We will next argue that for an arbitrary multimodc 
pure Gaussian state \^)abc the maximum E L ^ G between 
A and B is obtained if each party Cj performs homodyne 
detection of some quadrature xc^Oj ■ The proof is based 
on the reduction argument. Consider the multipartite 
state illustrated schematically in Fig. [5] Suppose that 
all parties Ck except for Cj perform their local projec- 
tions on pure Gaussian states described by POVMs |T]). 
This prepares the three-mode system ABCj in a pure 
Gaussian state with a fixed covariance matrix that does 
not depend on the measurement outcomes a,k of Ck, see 
Fig. [2 For a three- mode pure Gaussian state we proved 




=> 




FIG. 2: Reduction of iV-mode pure Gaussian state to three- 
mode pure Gaussian state by local projections on modes Ck- 



that El q is maximized if Cj makes a balanced homo- 
dyne detection of appropriately chosen quadrature. We 
can thus see that, irrespective of the measurements car- 
ried out by Cfe, the optimal choice for Cj is balanced 
homodyning. Note that Cj even does not need to learn 
the measurement outcomes ak of the other parties Ck^j- 
It suffices if A and B receive from all parties C/ the data 
ai . By local displacement operations they can then deter- 
ministically compensate for the resulting displacements 
of their modes which are linear in ai . The above argu- 
ment can be applied to any Ck which proves that the 
optimal measurement strategy must consist of balanced 
homodyning on each mode Ck- 

We now consider general mixed symmetric Gaussian 
states Q which are invariant under arbitrary permu- 
tation of modes. This implies that the covariance matrix 
has a highly symmetric form, 



7sym 



\ 



\ 



(13) 



where (3 and e denote symmetric 2x2 matrices. By 
means of local canonical transformations it is possible to 
simultaneously diagonalize both (3 and e fll| so without 
loss of any generality we may assume that (3 = diag(6, b) 
and e = diag(ei, £2)- We choose the logarithmic negativ- 
ity Ej^f as the measure of entanglement. We prove that 
maximum entanglement between two modes (labeled A 
and B) can be localized by homodyne detection of either 
x or p on all remaining N — 2 modes Cj . Consider the 
bipartite 2 x (N ~ 2) splitting AB\C. By means of uni- 
tary Gaussian transformation on modes Cj which can be 
physically realized by interference on an array of N — 3 
unbalanced beam splitters all modes Ck , k ^ 1, can be 
decoupled from A, B and C\ thereby effectively reducing 
the problem to three- mode case (Tl| . Similarly, interfer- 
ence of A with Bona balanced beam splitter decouples 
B from A and C\. In this way we obtain the equivalent 
representation shown in Fig.[T]where Uab represents mix- 
ing on a balanced beam splitter, mode B is initially in a 
mixed state with 75, i n = — e while the initial CM of 
modes A and C\ reads 



1AC U 



^2(N-2)e (3+(N-3)e 



(14) 
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FIG. 3: (a) Scheme for preparation of a three-mode Gaus- 
sian state, (b) Localizable entanglement El,g (solid curve) 
and El,ng (dashed curve) are plotted versus the squeezing 
parameter A. 



Projection of C\ onto state with CM 7c prepares mode 
A in state with CM ^ A ; m = (3 + e - 2(N - 2)e[f3 + 
(N — 3)e + 7c] _1 £- The logarithmic negativity Ej\f — 
max(0, — log 2 /i) where /1 is the minimum symplectic 
eigenvalue of CM of partially transposed state of modes 
A and B [16]. After some algebra one finds that [i 2 = 
min[eig(7 J 4 j i n i?7s j i n i? T )] where eig(^4) denotes eigenval- 
ues of a matrix A. In order to maximize Ej^f we have to 
minimize /i over all admissible jg - This could be done 
analytically and one can prove [151 ] that it is optimal to 
measure either the x or p quadrature of C\ , depending on 
the relation between e± and €2- This optimal measure- 
ment is a joint measurement on the original modes Cj 
that can be performed locally by measuring either x or 
p quadrature of each mode and then properly averaging 
the results, /i 2 is particularly simple for N = 3 when it 
reads /1 2 = (b - ei)(b - e 2 )(l + 2min(ei, e 2 )/b). 

Finally, we show that if we allow also non-Gaussian 
measurements on modes C then we can localize a higher 
amount of entanglement between modes A and B than 
one can localize by Gaussian measurements. This is 



yet another example of the well known fact that 11011- 
Gaussia n op erations are optimal for certain tasks such as 
cloning [TJ or partial estimation of coherent states [l8[ . 
We demonstrate the superiority of non-Gaussian mea- 
surements on a particular illustrative example. Consider 
the state preparation scheme shown in Fig. [3Ja). The 
three-mode state is generated from a two-mode squeezed 
vacuum state \\)bc m modes B and C by mixing mode 
Bona balanced beam splitter BS with mode A which is 
initially in a vacuum state. 

The optimal Gaussian measurement on C which maxi- 
mizes entanglement between A and B is a homodyne de- 
tection and the Gaussian localizable entanglement -El.g 

can be evaluated from formula ([6]) with tia — \ \J 1 \i — 

Suppose now that we would measure the number of 
photons n in mode C . With probability p n = (1 — A 2 ) A 2 ™ 
we would prepare in mode B n-photon Fock state 
which then impinges on a balanced beam splitter BS, c.f. 
Fig.[3Ja). The resulting state of A and B expre ssed in the 

Fock state basis reads, \ip n )AB = 5^72 YJl =0 \/GDl fc ' n_ 
k) ab- The entropy of entanglement of \ipn)AB is given by 



Sn 



2n Z 4 

k=0 







fnY 


(fc) log2 


— 

2" 





(15) 



and has to be evaluated numerically. The average entan- 
glement between A and B is then El,ng — J2^LoPnS n - 
Both El,g and El^ng are plotted in Fig.[3]Jb) as a func- 
tion of the squeezing parameter A. We can see that 
El,ng > q so the non-Gaussian measurement strat- 
egy outperforms the best Gaussian one. 
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